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Mean  Residual  Life:  Theory  and  Applications 

by 

Frank  Guess  and  Frank  Preschan 


1. 


Introduction  and  Summary. 

The  mean  residual  life  (MRL)  has  been  used  as  far  back  as  the  third  century 
A.D.  (cf.  Deevey  (1947)  and  Chiang  (1968)}.  In  the  last  two  decades,  however, 
reliabilists,  statistician*,  and  others  have  shown  intensified  interest  in  the  MRL 
and  derived  many  useful  results  concerning  it.  Given  that  a  unit  is  of  age  t,  the 
resaining  life  after  tiae  t  is  randon.  The  expected  value  of  this  random  residual 
life  is  called  the  mean  residual  life  at  time  t.  Since. the  MRL  is  defined  for  each 
tine  t,  we  also  speak  of  the  MRL  function.  (See  Section  2  for  a  sore  formal  defi¬ 
nition.) 

The  MRL  function  is  like  the  density  function,  the  moment  generating  function, 
or  the  characteristic  function:  for  a  distribution  with  a  finite  mean,  the  MRL 
completely  determines  the  distribution  via  an  inversion  formula  (e.g.,  see  Cox 
(1962),  Kotz  and  Shanbhag  (1980),  and  Hall  and  Wellner  (1981)).  Hall  and  Wellner 
(1981)  and  Bhattachar j ee  (1982)  derive  necessary  and  sufficient  conditions  for  an 
arbitrary  function  to  be  a  MRL  function.  These  authors  recommend  the  use  of  the 
MRL  as  a  helpful  tool  in  model  building. 

Not  only  is  the  MRL  used  for  parametric  modeling  but  also  for  nonparametric 
modeling.  Hall  and  Wellner  (1981)  discuss  parametric  uses  of  the  MRL.  Large  non- 
parametric  classes  of  life  distributions  such  as  decreasing  mean  residual  life 
(DMRL)  and  new  better  than  used  in  expectation  (NBUE)  have  been  defined  using  MRL. 
Barlow,  Marshall,  and  Pros chan  (1963)  note  that  the  DMIL  class  is  a  natural  one  in 
reliability.  Brown  (1983)  studies  the  problem  of  approximating  increasing  mean 
residual  life  (IMRL)  distributions  by  exponential  distributions.  He  mentions  that 
certain  IMRL  distributions,  "...  arise  naturally  in  a  class  of  first  passage  time 
distributions  for  Markov  processes,  as  first  illuminated  by  Keilson."  See  Barlow 
and  Proschan  (196S)  and  Hollander  and  Proschan  (1984)  for  further  comments  on  the 
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ncnparaaetric  use  of  >*tL. 

A  fascinating  aspect  about  WL  is  its  tremendous  range  of  applications.  For 
example,  Matson  and  Mells  (1961)  use  MRL  in  studying  bum-in.  Kuo  (1984)  presents 
further  references  on  MRL  and  bum- in  in  his  Appendix  1,  as  well  as  a  brief  history 
on  research  in  bum- in. 

Actuaries  apply  MRL  to  setting  rates  and  benefits  for  life  insurance.  In  the 
biomedical  setting  researchers  analyze  survivorship  studies  by  MUL.  See  Elandt- 
Johnson  and  Johnson  (1980)  and  Gross  and  Clark  (197S). 

Morrison  (1978)  mentions  IK®L  distributions  have  been  found  useful  as  models 
in  the  social  sciences  for  the  lifelengths  of  wars  and  strikes.  Bhattacharjee 
(1982)  observes  MRL  functions  occur  naturally  in  other  areas  such  as  optimal  dis¬ 
posal  of  an  asset,  renewal  theory,  dynamic  programming,  and  branching  processes. 

In  Section  2  we  define  more  formally  the  MRL  function  ind  survey  seme  of  the 
key  theory.  In  Section  3  we  discuss  further  its  wide  range  of  applications. 


2.  Theory  of  mean  residual  life. 


Let  F  be  a  life  distribution  (i.e.,  F(t)  *  0  for  t  <  0)  with  a  finite  first 
moment.  Let  F(t)  =  1  -  F(t).  X  is  the  random  life  With  distribution  F.  The  mean 
residual  life  function  is  defined  as 


(2.1) 


m(t)  *  E[X  -  t|X  >  t’j 
«  0 


for  F(t)  >  0 
for  F(t)  «  0 


for  ti  0.  Note  that  we  can  express  m(t)  «  /  d x*  J  du  when  F(t)  >0, 

0  F(t)  t  F(t) 

m 

If  F  also  has  a  density  f  we  can  write  m(t)  «  /  uf(u)du/F(t)  -t. 

t 
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Like  the  failure  rate  function  (recall  that  it  is  defined  as  r(t)  *f(t)/F{t) 
when  F(t)  >0),  the  MIL  function  is  a  conditional  concept.  Both  functions  are  con¬ 
ditioned  on  survival  to  tine  t. 

While  the  fa? lure  rate  function  at  t  provides  information  about  a  snail  inter¬ 
val  after  tine  t  ("just  after  t,"  see  p.  10  Barlow  and  Proschan  (1965)),  the  KRL 
function  at  t  considers  information  about  the  whole  interval  after  t  ("all  after 
t").  TTiis  intuition  explains,  the  difference  between  the  two. 

Note  that  it  is  possible  for  the  MRL  function  to  exist  but  for  the  failure 
rate  function  not  to  exist  (e.g.,  consider  the  standard  Cantor  ternary  function, 
see  Chung  (1974)  p.  12).  On  the  other  hand,  it  is  possible  for  the  failure  rate 
function  to  exist  but  the  MRL  function  not  to  exist  (e.g.,  consider  modifying  the 
Cauchy  density  to  yield  f(t)  «  ?  tjy  for  tiO).  Both  the  MRL  and  the  failure 

rate  functions  are  needed  in  theory  and  in  practice. 

When  a  and  r  both  exist  the  following  relationship  holds  between  the  two: 

(2.2)  a'(t)  «m(t)r(t)  -  1, 

for  a  differentiable.  See  Watson  and  Wells  (1961)  for  further  comments  on  (2.2) 
and  its  uses. 

If  the  failure  rate  is  a  constant  (>0)  the  distribution  is  an  exponential. 

If  the  MRL  is,  a  constant  (*0)  the  distribution  is  also  an  exponential. 

Let  ywE(X).  If  F(0)  »  0  then  o(0)  *■  u.  If  F(0)  >0  then  a(0)  «  y/F(0)  •  u.  For 
simplicity  in  discussions  and  definitions  in  this  section,  we  assume  F(0)  *0.  Let 
F  be  right  continuous  (not  necessarily  continuous) .  Knowledge  of  the  MRL  function 
completely  determines  the  reliability  function  as. follows: 
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where  F*1 (1)  d*{  sup{tjF(t)  <  1}. 

Cox  (1962)  assigns  as  an  exercise  the  demonstration  that  MRL  determines  the 
reliability.  Meilijson  (1972)  gives  an  elegant,  simple  proof  of  (2.3).  Kotx  and 
Shanbhag  (1980)  derive  a  generalized  inversion  formula  for  distributions  that  are 
not  necessarily  life  distributions.  Hall  and  Kellner  (1981)  have  an  excellent  dis¬ 
cussion  of  (2.3)  along  with  further  references. 

A  natural  question  to  ask  is:  what  functions  are  MRL  functions?  A  character¬ 
ization  is  possible  which  answers  this.  By  a  function  f  being  increasing  (decreas¬ 
ing)  we  mean  that  xiy  implies  f(x)  s£(i)  f(y). 


Theorem  2.1.  Consider  the  following  conditions: 

i.  a:  [0,  •)  -*-(0,  -). 

ii.  a((li>0. 

iii.  ■  is  fight  continuous  (not  necessarily  continuous), 

iv.  d(t)  d»^a(t)  ♦  t  is  increasing  on  [0,  •)• 

v.  When  there  exists  t_  such  that  a(tl)  lia  a(t)  »0,  then  a(t)  »0  holds 

t-t‘ 

for  t  c  [ tg ,  •)  .  Otherwise,  when  there  does  not  exist  such  a  tQ  with 

n(tl)  «  0 ,  then  /  -r  du ■  •  holds . 

0  0  “(u) 

A  function  n  satisfies  i  -  v  if  and  only  if  a  is  the  MRL  function  of  a  nondegenerate 
at  0  life  distribution. 

See  Hall  and  Kellner  .1981)  for  a  proof.  See  Bhattacharjee  (1982)  for  another 
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characterization.  Note  that  condition  ii  rules  out  the  degenerate  at  0  distribu¬ 
tion.  For  iv  note  that  d(t)  is  siaply  the  expected  time  of  death  (failure)  given 
that  a  unit  has  survived  to  tine  t.  Theorem  2.1  delineates  which  functions  can 
serve  as  MRL  functions,  and  hence,  provides  aodels  for  lifelengths. 

We  restate  several  bounds  involving  MRL  from  Hall  and  Wellner  (1981).  Recall 
a*  »  a  if  a  SO,  otherwise  a*«0. 

Theoren  2.2.  Let  F  be  nondegenerate.  Let  ur*EXrs»  for  r>l. 

i.  a(t)  £  (F'^l)  -  t)*  for  all  t.  Equality  holds  if  and  only  if 
F(t)  «F((F-|(1))‘)  or  1. 

ii.  a(t)  s  (y/F(t) )  -  t  for  all  t.  Equality  holds  if  and  only  if  F(t)  *0. 

iii.  a(t)  <  (u^F(t))l/r-  t  for  all  t. 

iv.  a(t)  2  (u  -  t)*/F(t)  for  t<F“1(l).  Equality  holds  if  and  only  if  F(t)*0. 

v.  a(t)  >[u-F(t)(wr/F(t))1/r]/F(t) -t  for  t<F“l(l). 

vi.  a(t)  s  (u  - 1)*  for  all  t.  Equality  holds  if  and  only  if  F(t)  *0  or  1. 

Various  nonparaaetric  classes  of  life  distributions  have  been  defined  using 
MRL.  (Recall,  for  simplicity  we  assuae  F(0)  *0  and  the  mean  is  finite  for  these  , 
definitions.) 

Definition  2.3.  PMRL.  A  life  distribution  F  has  decreasing  mean  residual  life  if 
its  MRL  a  is  a  decreasing  function. 

Definition  2.4.  NBUB.  A  life  distribution  F  is  new  better  than  used  in  expectation 


if  m(0)  a  m(t)  for  all  t  st  0 
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Definition' 2. b.  IDMRL.  A  life  distribution  F  has  increasing  then  decreasing  aean 
residual  life  if  there  exist  t  i  0  such  that  a  is  increasing  on  [0,  r)  and  decreas¬ 
ing  on  [r,  •). 

Each  of  these  classes  above  has  an  obvious  dual  class  associated  with  it*  i.e., 
increasing  aean  residual  life,  new  worse  than  used  in  expectation  (NWUE),  and 
decreasing  then  increasing  mean  residual  life  (D IMRL) ,  respectively. 

The  DMRL  class  models  aging  that  is  adverse  (e.g.,  wearing  occurs).  Barlow, 
Marshall,  and  Proschan  (lyt>3)  note  that  the  DMRL  class  is  a  natural  one  in  relia¬ 
bility.  See  also  Barlow  and  Proschan  (1965).  The  older  a  DMRL  unit  is,  the  shorter 
is  the  remaining  life  on  the  average.  Chen,  Hollander  and  Langberg  (1983)  contains 
an  excellent  discussion  of  the  vises  of  the  DMRL'  class. 

Bum- in  procedures  aie  needed  for  units  with  IMRL.  E.g.,  integrated  circuits 
have  been  observed  empirically  to  have  decreasing  failure  rates;  and  thus  they 
satisfy  the  less  restrictive  condition  of  IMRL.  Investigating  job  mobility,  social 
scientists  refer  to  IMRL  as  inertia.  See  Morrison  {1978)  for  example.  Brown  (1983) 
studies  approximating  IMRL  distributions  by  exponentials.  He  comments  that  certain 
IMRL  distributions,  "...  arise  naturally  in  a  class  of  first  passage  time  distribu¬ 
tions  for  Markov  processes,  as  first  illuminated  by  Keilson." 

Note  that  DMRL  implies  NBUE.  The  NBUE  class  is  a  broader  and  less  restrictive 
class.  Hall  and  Wellner  (1981)  show  for  NBUE  distributions  that  the  coefficient  of 
variation  cr/u  s  1,  where  a2  »  Var(X) .  They  also  comment  on  the  use  of  NBUE  in 
renewal  theory.  Bhattacharjee  (1984b)  discusses  a  new  notion,  age-smoothness,,  and 
its  relation  to  NBUE  for  choosing  lifa  distribution  models  for  equipment  subject  to 
eventual  wear.  Note  that  burn-in  is  appropriate  for  NWUE  units. 

For  relationships  of  DMRL,  IMRL,  NBUE,  and  NWUE  with  other  classes  used  in 


-  7  - 


reliability  see  the  survey  paper  Hollander  and  Proschan  (1964). 

The  IDMRL  class  models  aging  that  is  initially  beneficial,  then  adverse.  Sit¬ 
uations  where  it  is  reasonable  to  postulate  an  IDMRL  model  include: 

i  Length  of  time  employees  stay  with  certain  coopanies:  An 
employee  with  a  company  for  hnir  years  has  more  tJae  and  career 
invested  in  the  company  than  an  employee  of  only  two  months. 

The  MRL  of  the  four-year  employee  is  likely  to  be  longer  than 
the  MRL  of  ths  two-month  eaplovee.  After  this  initial  IMRL 
(this  is  called  "inertia"  by  social  scientists),  the  processes 
of  aging  and  retirement  y*eld  a  CMRL  period. 

ii  Life  lengths  of  humans;  High  infant  mortality  explains  the 
initial  IMRL.  Deterioration  and  aging  explain  the  later  DMRL 
stage. 

See  Guess  (1984)  and  Guess,  Hollander,  and  Pros chan  (1983)  for  further  examples  and 
discussion.  Bhattacharjee  (1983)  comments  that  Gertsbakh  and  Kordonskiy  (1969) 
graph  the  MRL  function  of  a  lognormal  distribution  that  has  a  "bath-tub"  shaped  MRL 
(i.e.,  DIMRL). 

Hall  and  We liner  (1981)  characterize  distributions  with  MRL's  that  have  linear 
segments.  They  use  this  characterization  as  a  tool  for  choosing  parametric  models. 
Morrison  (1978)  investigates  linearly  IMRL.  He  states  and  proves  that  if  F  is  a 
mixture  of  exponential  then  F  has  linearly  IMRL  if  and  only  if  the  mixing  distribu¬ 
tion,  say  G,  is  a  gamma.  Howell  (1984)  studies  and  lists  other  references  on  lin¬ 
early  DMRL. 

In  renewal  theory  MRL  arises  naturally  also.  For  a  renewal  process  with  un<ier- 

cm, 

lying  > istribution  F,  let  G(t)  •  (jF(u)du)/u.  G  is  the  limiting  distribution  of  both 

t 
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the  forward  and  the  backward  recurrence  tines.  ee  Cox  (l'*62)  for  sore  details. 

Also  if  the  renewal  process  is  in  equilibrium  then  G  is  the  erect  distribution  of 

the  recurrence  tines.  G(t)  *  (m(t)F(t))/y.  The  failure  rate  of  G,  r^,  is  inversely 

related  to  the  MRL  of  F,  Op.  ,  I.e.,  r^(t)  *  l/r»p(t).  Note,  however,  that 

rp(t)  * l/aptt)  is  usually  the  case.  See  Hall  and  Kellner  (1S81),  Rolski  (1975), 

Meilijson  (1972),  and  ttatson  and  Hells  (1961)  for  related  discussions. 

Kot2  and  Shanbhag  (1980)  establish  a  stability  result  concerning  convergence 

of  an  arbitrary  sequence  of  MRL  functions  to  a  limiting  MRL  function.  (See  also' 

Bhattacharjee  (1982)).  They  show  an  analogous  stability  result  for  hazard  aeasures. 

(When  the  failure  rate  for  F  exists  and  Vp  is  F’s  hazard  measure,  then  vp(B) 

*  J  r_(t)dt  for  B  a  Borel  set.)  Their  results  imply  that  MkL  functions  can  provide 
B  ■ 

more  stable  and  reliable  information  than  hazard  measures  when  assessing  noncontin- 
uous  distributions  from  data. 

In  a  multivariate  setting,  Lee  (1985)  shows  the  effect  of  dependence  by  total 
pdsitivity  on  MRL  functions. 


3.  Applications  of  mean  residual  life. 


A  mean  is  easy  to  calculate  and  explain  to  a  person  not  necessarily  skilled 
in  statistics.  To  calculate  the  empirical  MRL  function,  one  does  not  ne**i  calculus. 
Details  of  computing  the  empirical  MRL  follow. 

Let  X^,  •  ••*  xn  be  a  random  sample  from  F.  For  simpler  initial  notation, 

we  assume  first  no  ties.  Later  we  allow  for  ties.  Order  the  observations  as 


(3.1) 


Xln<*2n<  *"*  <Xnri* 


Let  X^‘0.  The  empirical  MRL  function  is  defined  as 
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l  <*in-s) 

i  *k+I  ln 


for  t£  fXkn’  X(k+l)n)l 


and  k  «  C ,  1 . n  -  1 .  a  (t)*0  for  t  2  X  . 

nv  *  nn 

Note  that  (3.2)  is  simply 


(3.3) 


.  .  Total  time  on  test  observed  after  t 
an^  *  Number  of  units  observed  after  t 


.  *  xi 

The  empirical  KRL  function  at  0,  m  (0)  »X  *  (  £  X.)/n,  is  just  the  usual  sample 

n  n  •  •  x 

1*1 

mean  when  no  unit  fails  at  time  0.  If  a  unit  fails  at  0  then  m  (0)  >X  . 

n  n 

If  ties  jxist  let 


(3.4) 


be  the  distinct  ordered  times  of  failure. 


(3.S)  n.  « number  of  observed  failures  at  time  X.»,  s.  =  n-  £  n. 

1  j  *0  3 

for  i*0,  1,  t<  n.  Note  that  n.<0,  i*l,  ...,  t,  while  n^'O  is  allowed. 


(3.6) 


®n(t)  *X  ^  ls  for  te  X(k+lK5' 


■n(t)  *  0 


for  taX. 


for  k*0r  1,  ....  4*1.  Note  that  (3.6)  is  simply  notation  for  (3.3). 
We  illustrate  in  the  following  example, 


3.1.  Bjerkedal  (1960)  studies  the  lifelengths  of  guinea  pigs  injected  with 
different  amounts  of  tubercle  bacilli.  Guinea  pigs  are  known  to  have  a  high  sus¬ 
ceptibility  to  human  tuberculosis,  which  is  one  reason  for  choosing  this  species. 


We  describe  the  only  study  (M)  in  which  animals  in  a  single  cage  are  under  the  same 
regimen.  The  regimen  number  is  the  common  log  of  the  number  of  bacillary  units  in 
.5  ml  of  the  challenge  solution,  e.g.,  regimen  4.3  corresponds  to  2.2 x  1G4  bacil¬ 
lary  units  per  .5  ml  (log1Q(2.2 * 104)  *4.342).  Table  3.1  presents  the  data  from 
regimen  5.5  and  the  empirical  MRL. 

Graphs  of  MRL  provide  useful  information  not  only  for  data  analysis  but  also 

for  presentations.  Commenting  on  fatigue  longevity  and  on  preventive  maintenance, 

» 

Gertsbakh  and  Kordonskiy  (1969)  recommend  the  MRL  function  as  another  helpful  tool 
in  such  analyses.  They  graph  the  MRL  for  different  distributions  (e.g.,  Weibull, 
lognormal,  and  gamma) .  Hall  and  Wellner  (1979)  graph  the  empirical  MRL  for 
Bjerkedal's  (1960)  regimen  4.3  and  regimen  6.6  data.  Bryson  and  Siddiqui  (1969) 
illustrate  the  graphical  use  of  the  empirical  MRL  on  survival  data  from  chronic 
granulocytic  leukemia  patients.  Using  the  standard  Kaplan-Meier  estimator  (e.g., 
see  Lawless  (1982),  Nelson  (1982),  or  Miller  (1980)),  Chen,  Hollander,  and  Langberg 
(1983)  graph  the  empirical  MRL  analogue  for  censored  lifetime  data. 

Gertsbakh  and  Kordonskiy  (1969)  note  that  estimation  of  MRL  is  more  stable 
than  estimation  of  the  failure  rate.  Statistical  properties  of  estimated  means  are 
better  than  those  of  estimated  derivatives  (which  enter  into  failure  rates). 

Yang  (1978)  shows  that  the  empirical  MRL  is  uniformly  strongly  consistent. 

She  establishes  that  suitably  standardized,  converges  weakly  to  a  Gaussian  pro¬ 
cess.  Hall  and  Wellner  (1979)  require  less  restrictive  conditions  to  apply  these 
results.  They  derive  and  illustrate  the  use  of  simultaneous  confidence  bands  for 
m.  Yang  (1978)  comments  that  for  t>0,  mn(t)  is  a  slightly  biased  estimator. 

Specifically,  E(mn(t))  * m(t)  (1  -  Fn(t)) .  Note,  however,  that  liaE(m  (t))  »m(t) . 

" 

Thus,  for  larger  samples  »n(t)  is  practically  unbiased.  See  also  Gertsbakh  and 
Kordonskiy  (19*9). 


TABLE  3.1 


Empirical  mean  residual  life  in  days  at  the  unique 
times  of  death  for  the  72  guinea  pigs  under  regimen  5.S. 
(We  include  the  empirical  MRL  at  time  0  also.) 


Number  of  Ties  Time  of  Death  Empirical  MRL 


n. 

X. 

m  (X.  : 

i 

m 

nv  m 

0 

0 

141.85 

1 

43 

100.24 

1 

45 

99.64 

1 

53 

92.97 

2 

56 

92.66 

1 

57 

93.05 

1 

58 

93.46 

1 

66 

86.80 

1 

67 

87.16 

1 

73 

82.47 

1 

74 

82.80 

1 

79 

79.10 

2 

80 

80.79 

3 

81 

84.15 

1 

82 

84.69 

2  , 

83 

86.90 

1 

84 

87.59 

1 

88 

85.26 

1 

89 

85.98 

2 

91 

87.55 

2 

92 

90.40 

1 

97 

87.34 

2 

99 

89.40 

2 

100  • 

92.83 

1 

101 

94.18 

3 

102 

100.94 

1 

103 

102.80 

1 

104 

104.79 

1 

107 

104.88 

1 

108 

107.13 

1 

109 

109.55 

1 

113 

109.07 

1 

114 

111.79 

1 

118 

111.64 

1 

121 

112.67 

-  . _  ’  _ '  _ ; _  •  — _ ; _ f  'i 
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Number  of  Ties 


Table  3.1  (Continued) 


Tine  of  Death 


X. 

m 

n  (X.  ) 
nv  m 

123 

114.92 

126 

116.40 

128 

119.17 

137 

114,96 

138 

119.14 

139 

123.76 

144 

124.70 

145 

130.21 

147 

135.33 

156 

133.76 

162 

135. 75 

174 

132.03 

178 

137.14 

179 

146.62 

184 

153.42 

191 

159.73 

198 

168.00 

211 

172.22 

214 

190.38 

243 

184.43 

249 

208.17 

329 

153.80 

380 

128.50 

403 

140.67 

511 

49.00 

522 

76.00 

598 

0.00 
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Yang  (1977)  studies  estimation  of  the  HRL  function  when  the  data  are  randomly 
censored.  Joe  and  Proschan  (1981)  develop  isotonized  estimators  of  MRL  and  of  the 
life  distribution  under  the  assumption  of  DMRL.  Both  the  complete  data  model  and 
the  randomly  censored  data  model  ere  treated. 

For  parametric  modeling  Hall  and  Wellner  (1981)  use  the  empirical  MRL  plot. 
They  observe  that  the  empirical  MR&  function  is  a  helpful  addition  to  other  life 
data  techniques,  such  as  total  time  on  test  plots,  empirical  (cumulative)  failure 
rate  functions,  etc.  The  MRL  plot  detects  certain  aspects  of  the  distribution  more 
readily  than  other  techniques.  See  Hall  and  Wellner  (1981),  Hall  and  Wellner 
(1979),  and  Gertsbakh  and  Kordonskiy  (1989)  for  further  comments. 

When  a  parametric  approach  seems  inadvisable,"  the  MRL  function  can  still  be 
used  as  a  nonparametric . tool .  Broad  classes  defined  in  terms  of  MRL  allow  a  more 
flexible  approach  while  still  incorporating  preliminary  information.  For  example, 
to  describe  a  wear  process,  a  DMRL  is  appropriate.  When  newly  developed  components 
are  initially  produced,  many  may  fail  early  (such  early  failure  is  called  infant 
mortality  and  this  early  stage  is  called  the  debugging  stage).  Another  subgroup 
tends  to  last  longer.  Depending  on  information  about  this  latter  subgroup,  we 
suggest  IMRL  (e.g.,  lifelengths  of  integrated  circuits)  or  IDMRL  (e.g.,  more  com¬ 
plicated  systems  where  there  are  i.ifant  mortality,  useful  life,  and  wear  out 
stages). 

Objective  tests  exist  for  these  and  other  classes  defined  in  terms  of  MRL. 
E.g.,  see  Hollander  and  Proschan  (1984)  and  Guess,  Hollander,  and  Proschan  (1983). 

To  describe  "bum- in”  the  MRL  is  a  natural  function  te  use.  Kuo’s  (1984) 
Appendix  1  presents  an  excellent  brief  introduction  to  bum- in  problems  and  appli¬ 
cations  of  MRL. 

Actuaries  apply  MRL  to  setting  rites  and  benefits  for  life  insurance.  In  the 
biomedical  setting  researchers  analyze  survivorship  studies  by  MRL.  For  example, 
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see  Elandt-Johnson  and  Johnson  (1981)  and  Gross  and  Clark  (1975). 

Social  scientists  use  IMRL  for  studies  on  job  mobility,  length  of  wars,  dura¬ 
tion  of  strikes,  etc.  See  Morrison  (1978). 

In  economics  MRL  arises  also.  Bhattacharjee  and  Krishnaji  (1981)  present 
applications  of  MRL  for  investigating  landholding.  Bhattacharjee  (1984a)  uses  N8UE 
for  developing  optimal  inventory  policies  for  perishable  items  with  random  shelf 
life  and  variable  supply. 

Bhattacharjee  (1982)  observes  MRL  functions  occur  naturally  in  other  areas 
such  as  optimal  disposal  of  an  asset,  renewal  theory,  dynamic  programming,  and 
branching  processes. 
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